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Occurrence of strong or mass synchronization of a large number of neuronal populations in the 
brain characterizes its pathological states. In order to establish an understanding of the mechanism 
underlying such pathological synchronization we present a model of coupled populations of phase 
oscillators representing the interacting neuronal populations. Through numerical analysis, we discuss 
the occurrence of mass synchronization in the model, where a source population which gets strongly 
synchronized drives the target populations onto mass synchronization. We hypothesize and identify 
a possible cause for the occurrence of such a synchronization, which is so far unknown: Pathological 
synchronization is caused not just because of the increase in the strength of coupling between the 
populations but also because of the strength of the strong synchronization of the drive population. 
We propose a demand-controlled method to control this pathological synchronization by providing a 
delayed feedback where the strength and frequency of the synchronization determines the strength 
and the time delay of the feedback. We provide an analytical explanation for the occurrence of 
pathological synchronization and its control in the thermodynamic limit. 

PACS numbers: 



Strong synchronization of several groups of neu- 
ronal oscillators occurs during pathological states 
like Parkinson's tremor, epileptic seizures, etc. 
This synchronization first occurs locally in a par- 
ticular region of the brain due to synchroniza- 
tion of a group of neurons. This group then 
acts as a pacemaker triggering synchronization to 
various other neuronal populations resulting in a 
mass synchronization. Deep brain stimulation is 
a standard therapy for patients with this disor- 
der. However since the actual cause and mech- 
anism of the occurrence of mass synchronization 
in the brain is not yet fully known, there is al- 
ways a need to understand the same in order to 
be able to control it effectively. In this article 
we attempt to gain this understanding by study- 
ing a system of coupled populations of phase os- 
cillators representing interacting neuronal popu- 
lations in the brain. We find the occurrence of 
mass synchronization in such a system where the 
occurrence of synchronization in the source popu- 
lation triggers synchronization in the target pop- 
ulations also. We hypothesize a possible cause 
for the occurrence of mass synchronization in the 
brain, which is so far not known. We propose a 
demand controlled method to control this mass 
synchronization by providing a delayed feedback. 
We provide numerical evidence and analytical ex- 
planation for the occurrence and control of neu- 
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ronal mass synchronization. 



I. INTRODUCTION 

Synchronization is an ubiquitous phenomenon and a 
topic of active research in recent years [ll-Q . Given that 
it is not always a desirable phenomenon, synchrony con- 
trol mechanism are always required. For example, in the 
brain, when synchronization supports cognition via tem- 
poral coding of information [4, 5] , it is desirable while at 
the same time it is undesirable when synchronization of 
a mass of neuronal oscillators occurs at a particular fre- 
quency band resulting in pathologies like trauma, Parkin- 
son's tremor and so on. Other examples include lasers 
and Josephson junction arrays @, emission of microwave 
frequencies by coupled spin torque nano oscillators 0] 
where synchronization is desirable, while in the cases of 
epileptic seizures [8i], Parkinson's tremor 0], event re- 
lated desynchronization [l^, UII, or pedestrians on the 
Millennium Bridge it is undesirable. 

Pathological strong (or mass) synchronization occurs 
when several groups of neurons enter into synchroniza- 
tion and the strength of the synchronization is large 
because of the large number of neuronal networks that 
are engaged in synchrony. In particular, tremors in the 
case of Parkinson's disease (PD), essential tremor and 
several other neurological diseases like epileptic seizures 
are caused by mass synchronization of oscillating neu- 
ronal networks. Synchronization in such cases may also 
be called abnormal synchronization, since under normal 
conditions mass synchronization does not occur. Clus- 
ters of neurons fire in a synchronized manner with a fre- 
quency similar to that of the tremor and act like a pace- 
maker driving the pre-motor and motor cortical areas of 
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the brain causing mass synchronization |13l . Il4| . In fact 
it has been shown that the local field potentials in the 
thalamic and the sub-thalamic nucleus (STN) drive the 
tremor in the limb In the case of healthy subjects 

these neurons in the STN oscillate in a desynchronized 
manner [l8l |. Deep brain stimulation is a standard ther- 
apy for such patients with such tremors where a train 
of strong high frequency electric pulse is administered 
to the areas like STN where there is mass synchroniza- 
tion. In spite of the great amount of interest and research 
focussed on pathological tremors, the actual mechanism 
that cause these tremors are still unknown [T9l - [23j . Fur- 
ther, controlling of pathological tremors is also a problem 
which has recently gained much attention [23t22| ■ 

Given that the interaction mechanism and the topolog- 
ical structure of the neuronal networks are highly com- 
plex and are not completely known, and that there is a 
lack of technology to measure synchronization in every re- 
gion of neuronal network, it is highly complicated to gain 
a complete understanding of the underlying mechanism 
in these pathological states. In fact, the actual mecha- 
nism for synchronous oscillations in Parkinson's disease 
remains unknown. Nevertheless, one can always repre- 
sent the system of interacting neurons as a mathematical 
model and simulate the challenging case of a pathological 
synchronization to understand the underlying dynamical 
aspects. With this motive, in this paper, we hypothesize 
the possible cause for the occurrence of mass synchro- 
nization in PD tremors and also provide a method to 
control it. 

The plan of the paper is as follows: In the following 
Sec. In] we introduce a model of coupled populations 
of phase oscillators which represents interacting popu- 
lations of neurons. We present our numerical finding of 
the occurrence of mass synchronization where the source 
population drives the target in Sec. IIIII In Sec. ITVl we 
propose a method to control the occurrence of pathologi- 
cal synchronization by providing a delayed feedback. We 
analytically explain the occurrence of mass synchroniza- 
tion and the method to control it in Sec. |Vl Finally in 
Sec. IVII we present our conclusions. 

II. THE MODEL 

We consider a model of interacting populations of neu- 
rons which is represented by populations of coupled phase 
oscillators. One of the populations is considered to be a 
source which acts as a pacemaker and possibly drives the 
other populations (named as targets) onto mass synchro- 
nization. The model equations read as 

N' N^'=N 
»7' = 1 j=l 

i = l,2,...,N^, =N,{1) 

where i] = 1,2,... TV', N' is the number of populations 
and Nri' is the number of neurons in population 77' which 



is assumed to be the same in all the populations, that 
is equal to N. a^^i is the strength of the coupling be- 
tween the neurons in r/ and those in r/. Note that 
r] — rj' (1,2, ... N') corresponds to interconnection among 

(v) 

the neurons of the same population. Here is the nat- 
ural frequency of the kth neuron in the population ij and 
< \oirjr]'\ < is the phase lag. In reality every oscil- 
lator has different a, but we have taken them to be the 
same for all of them, for simplicity. 

For a better understanding, a schematic diagram of 
the model is given in Fig. [T]for A^' = 2, where the source 

population is represented by a dotted circle {0^^'') and the 

(2) 

target is represented by a solid circle {9- ). In normal 
states it is the low frequency signal that drives the high 
frequency signal [H, [2^ , but in the case of pathological 
synchronization, high frequency signal from the source 
drives the target populations that are oscillating at (rel- 
atively) low frequencies on to synchronization. Therefore 
we have chosen ai2 — fJ-cr2i, where < fi < I, since the 
coupling from the source to the target (cr2i) is stronger 
than the coupling from the target to the source (0-12) in 
pathological states. Oscillators inside the source and the 
target populations are coupled among themselves with 
strengths an and (T22, respectively. In order to quantify 
synchronization within a population, we can define the 
order parameter as 

1 ^ 

z, = r,e^^<''>^l5:e-r. (2) 
i=i 

When r,, — 1, there is complete synchronization within 
the rj population since in this state the phases of all the 
oscillators in that population are the same. When 
takes values in between and 1 there is partial desyn- 
chronization and when r,, — there is complete desyn- 
chronization in the 77 population. The time average of r,, 
is given by 

1 

R,j =< r,j >= y y "^v^t, (3) 

which is used to characterize the occurrence of strong 
synchronization in a population. Numerically, strong 
synchronization is characterized by i?,, > 0.8 and desyn- 
chronization is characterized by Rj^ < 0.3 for T = 10^ 
time units. When strong synchronization occurs in all 
the populations due to the strong coupling strength in 
the source it is called mass synchronization. From equa- 
tions H]) and ^ the model equations can be rewritten 
as 

N' 

^/"^ = c^i") - J2 ^m'W ^HO^:"^ - V'^"'^ + ar,,'). (4) 

7,' = ! 

In order to find out the cause for the occurrence of mass 
synchronization and also to find ways to control it, we nu- 
merically simulate system Q using Runge-Kutta fourth 
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FIG. 1: Schematic representation of system ^ for A^' = 2. 
Sj-^' is the source population and S'^' is the target. The cou- 
pling strengths within the populations are quantified by the 
parameters an and a22- The coupling strengths from the 
source to the target and the target to the source are quanti- 
fied by the parameters ai2 = fj,a2i and (721 , respectively. 
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order routme with a time step of 0.01. We have set 
N — 1000 and have assumed a Lorentzian distribution 
for the oscillator frequencies given by 
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-7r, 
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where 7 is the half width at half maximum and w,, is the 
central frequency. We have discarded transients of the 
order of 10^ in our simulations. Simulations lasted for 
2 X 10^ steps and the results are shown for a small window 
of time (which we label to start with in the figures 
for convenience) towards the end of the total simulation 
time. The initial phases are considered to be randomly 
distributed between and 27r. In the following sections 
we will discuss how the occurrence of synchronization 
in the source population causes mass synchronization in 
the target populations and will also present a method to 
control this mass synchronization. 

The model we consider here is of mass model type and 
focusses on the macroscopic properties such as the syn- 
chronization of large population of neurons. The model 
therefore does not include a detailed microscopic descrip- 
tion to incorporate the effects of hierarchal connections, 
spiking activities of individual neurons and populations, 
the role of neurotransmitters and so on [30l - l34| . We make 
use of the macroscopic modeling approach to strike a bal- 
ance between abstract and detailed modeling while keep- 
ing the number of variables and parameters to a reason- 
able minimum. Thus we attempt to make a qualitative 
concurrence between the model results and the occur- 
rence of pathological mass synchronization. 



III. OCCURRENCE OF MASS 
SYNCHRONIZATION 

In order to illustrate the occurrence of mass synchro- 
nization, we have taken iV' = 2, so that we have a sys- 
tem of one target and one source populations, each hav- 
ing 1000 oscillators. We have verified our results for the 
A^' = 3 case also (see Fig. S] below) . To start with, 
we consider a state when both the source and the target 



FIG. 2: The time evolution of the phases in the source (panels 
(a) and (c)) and the target (panels (b) and (d)) populations 
are plotted for A'^ — 1000. Panels (a) and (b) represent the 
initial desynchronized state characterized by -Ri — 0.2 and 
R2 = 0.3 and panels (c) and (d) represent the occurrence of 
strong synchronization in the target {R2 — 0.98) when the 
source population s'^' gets synchronized (_Ri = 0.99). Here 

(722 = 0.1, (721 ~ 1.5, fl — 0.3, 71,2 ~ 0.05, CJl = 1.5, LjJ2 = 0.5, 

aij — 0, i,j — 1, 2, and (7ii — 0.1 in (a) and (b) and an — 2.5 
in (c) and (d). 



populations are desynchronized (cth — 022 — 0.1) with a 
very low value of the average order parameters R\ — 0.2 
and i?2 = 0.3, even with a strong coupling between the 
source and the target {021 = 1.5 and fx = 0.3). Upon 
increasing the coupling strength within the source popu- 
lation, that is (Til = 2.5, we find that strong synchroniza- 
tion occurs in the source (characterized by i?i = 0.99) 
which in turn induces strong synchronization in the tar- 
get (characterized by R2 — 0.98). Note that this strong 
synchronization occurs in the target due to its coupling 
with the source, that is (T21, and not because of the cou- 
pling within itself, since tT22 = 0.1 only. This can be seen 
from Fig. [2] where we have plotted the time evolution 
of the phases in the source (panels (a) and (c)) and the 
target (panels (b) and (d)) populations. Panels (a) and 
(b) show the initial desynchronization state and (c) and 
(d) depict the occurrence of mass synchronization due to 
the increase in the coupling (an) within the source pop- 
ulation. Fig. [3] shows the snapshots of the distribution of 
the phases in the (x^, ?;i)=(cos 0^, sin 0;) plane where the 
panels (a)-(d) correspond to those in Fig. [21 In panels 
(a) and (b) the oscillators in the source and the target 
populations are desynchronized initially and hence the 
phases are distributed throughout the unit circle. Pan- 
els (c) and (d) depict the distribution of phases in the 
source and the target populations, respectively, in the 
mass synchronization state; the phases occupy a rather 
narrow region on the unit circle. 

Thus it becomes evident that for the onset of patholog- 
ical synchronization the neurons in the source population 
are strongly coupled among themselves. Due to this, the 
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FIG. 3: (a)-(d) are the snapshots of the phase portraits on 
the {xi,yi) plane corresponding to Fig. [2](a)-(d), where Xi — 
cos 6i and yi = sin 6i . 
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FIG. 4: The time average of r,, for increasing cru is plotted 
for two cases, N' — 2 (dotted lines) and N' = 3 (solid lines). 
Here cr,,,, = 0.05, cr,,i = 1.5, air, = i-icrni, for rj = 2,3, = 
0.2, = 7r/2 - 0.1, i,j = 1,2,3, 71,2,3 = 0.05, uji = 1.5, 
'»J2,3 = 0.5, and a23 = 0-32 — 0.01. The regions DS and 
SS denote the desynchronization and strong synchronization 
states, respectively. 



occurrence of a strong synchronization in the source pop- 
ulation is propagated onto the target populations because 
of an increase in the order parameter ri of the source, 
see equation (|4]). As a consequence a2iri increases faster 
than /i(T2ir2 and hence the source drives the target and 
makes it to get synchronized with itself. Under normal 
conditions, however an does not attain higher values to 
cause strong synchronization in the source. Thus for the 
occurrence of mass synchronization not only the coupling 
strength between the populations is crucial but also the 
strength of the synchronization of the source (character- 
ized by ri or i?i) plays a vital role. This is evident from 
Fig. m where we have plotted the time average of (i?^) 
for varying coupling strength an of the source for the 
two cases, N' = 2 and TV' = 3. 

For the case iV' = 3 we have considered three pop- 
ulations of coupled phase oscillators each having 1000 
oscillators with cr^^ = 0.05, cr^i — 1.5, air, — /^Cjyi, for 

= 2,3, n ^ 0.2, 71,2,3 = 0.05, wi = 1.5, ^2,3 = 0.5, 
and (723 = 0-32 = 0.01. One of the populations is the 
source (with average order parameter Ri) and there are 
two target populations (with average order parameters 
i?2 and R3) in the system. For both the cases N' = 2 
and TV' = 3, as cth increases the strength of the syn- 
chronization of the source increases and as a consequence 
synchronization in the target populations also increases. 
Interestingly, as one can see from the figure, for a given 
an, the target population achieves strong synchroniza- 
tion with the source even before the source population 
achieves complete synchronization. Our model results 
prove to be valid for a system of large number of coupled 
populations also and we have presented here the case of 
three populations for illustration. 



IV. METHOD TO CONTROL MASS 
SYNCHRONIZATION 

In order to control the occurrence of mass neuronal 
synchronization we apply a delayed feedback force to the 
populations in which case the model equation ([1]) be- 
comes 

N' 

= ^1") _ ^ a,,,r,, sin(0('') - i^^^'^ + a,,,) 

r,' = l 

±FRrSm{9^'''' -(j)r), i^l,2,...,N, (6) 

where Rr — R{t — t) and (j)r — (t){t — r) and Z{t) = 

E^^Li^r,' = Re'^ = i^E^Li'^r^'e''^''' is the global or- 
der parameter. The feedback can be either positive or 
negative as denoted by the ± sign in the above equa- 
tion since the dynamics of the system does not change 
qualitatively upon using either -f or — signs [s^. Also, 
one can note that the ± sign can be interchanged upon 
replacing cf) with — tt. 

The global order parameter R{t) measures the strength 
of synchronization of all the populations (combined) in 
the system while the order parameter r,, measures the 
strength of synchronization of the 77th population. The 
schematic representation of this mass synchronization 
control set up is shown in Fig. [5l The strength of the 
global synchronization R{t) and the corresponding fre- 
quency {— d(j)/dt) are measured from the populations 
of neurons denoted by filled circles. A delayed feedback 
of strength F and a time delay of r is fed back to the 
neuronal populations in order to control mass synchro- 
nization. This method is effective since one can have 
control over certain parameters like time delay and the 
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FIG. 5: (Color online) Schematic representation of the con- 
trol set up. Filled circles represent the neuronal populations. 
The strength of the global synchronization (R) and the corre- 
sponding frequency (Q) are measured. A delayed feedback of 
strength F and a time delay of r is fed back to the neuronal 
populations in order to cause desynchronization. 




FIG. 6: (a) and (b) show the occurrence of desynchronization 
for r = 1.0 and F = 5.5 in the source (s'^') and the tar- 
get populations, respectively due to delayed feedback, 
following panels (c) and (d) of Fig. [2] Here N=256. 



strength of the feedback. There is also no need to know 
about the interactions among the neurons or the charac- 
teristics of individual neuronal oscillators, but all we need 
is the macroscopic information. Further, the feedback is 
demand controlled because F and r are fixed depending 
upon the values of R and fi. In addition, there is no 
need to seek for a new controlling signal that stimulates 
the neuronal populations, but the measured pathological 
signal is used to control the same. We choose to provide 
delayed feedback because it can affect the activity of the 
individual neurons by shifting their phases. 

The numerical simulation of equation ^ showing the 
control of mass synchronization is given in Fig. |6] which 
follows the state of mass synchronization shown in Fig. [2] 
(panels (c) and (d)) after applying the delayed feedback 
(r = 1.0 and F — 5.5). Due to the delayed feedback both 
the source (panel (a)) and the target (panel (b)) popula- 
tions are desynchronized. The corresponding distribution 
of phases on the {xi,yi) plane is shown in Fig. [7] where 
the phases are distributed throughout the unit circle de- 
noting the desynchronization states of the source (panel 
(a)) and the target (panel (b)) populations. For numer- 
ical simulations including delayed feedback we have set 

= 256. We have also checked that the results are size 
independent by calculating them for =128 and 512. 

In the following section we analytically show the oc- 
currence of mass synchronization and its control with 
delayed feedback in support of the above numerical find- 
ings as follows: (i) We will show that when the source is 
desynchronized the target population will also be desyn- 
chronized. (ii) Then when strong synchronization occurs 
in the source the same is induced in the target also, (iii) 
In the strong synchronization state, when delayed feed 
back is introduced the synchronization is broken for suit- 
able values of the strength of the feedback and the value 
of the time delay. 



V. ANALYTICAL EXPLANATION 

In order to analytically explain the occurrence of mass 
synchronization and its control, we analyze system ([T]) 
in the continuum limit N ~> oo. In this limit, a prob- 
ability density for the oscillator phases can be defined 
as p'^'^^ {uj, 6, t), which describes the number of oscillators 
with phases within [6, 9 + d9] and natural frequencies be- 
tween [Lij,u! + duj] at time t. This distribution p^'^\u!,9,t) 
obeys the evolution equation 

^ = -^(p''')«'''^^ = l,2,...,iV' (7) 

where from equation (|4]) (in the continuum limit) the 
velocity u^''' is given by 
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7? = 1,2,...,7V'. (8) 

Also, in the continuum limit the order parameter ^ 
for the population rj becomes 
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Now we rewrite the velocity in terms of the order param- 
eters to get 

N' , 
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FIG. 7; (a) and (b) are the snapshots of the phase portraits 
on the {x,y) plane corresponding to Fig. [5] (a) and (b) for 
N=256. 



Upon introducing the ansatz of Ott and Antonsen [36[ 
we consider a special class of the density functions in the 
form of Poisson kernel as 



^ L ^k=l 

,77 = 1,2,. ..,7V', 



(11) 



where giuj'^'^'^ is the distribution of the oscillator frequen- 
cies. Here c.c denotes the complex conjugate of the pre- 
ceding term. The main advantage of using this ansatz is 
that it has the same function a,, in all the Fourier har- 
monics, except that a,, is raised to the nth power in the 
nth harmonic. Thus one can reduce the infinite set of 
dynamical equations exactly into a set of finite number 
of equations. Using the above expression ([TT|) for p^^^ 
into the evolution equation ([7]), we get the equation for 
Ctrl as 



N' 



r)' = l ^ 



e '"""'z^/a^ I, 



Now we express the amplitude equation (1141) in polar 
coordinates = rf^e^i as (see equation ^) 



1 — r: 



2 N' 



r-r, = -Jv^rj + { — ^) E cr,,,,.r^/ cos(V',, - V^^' + a,,r,')> 



1 + rf ^' 

77= l,2,...,iV'. (15) 

Using these equations we will explain the occurrence of 
pathological synchronization when the source is synchro- 
nized, and its control by applying delayed feedback, in 
the following subsections. 

A. Occurrence of mass synchronization 

As noted earlier in Sec. Ill, we have numerically found 
the occurrence of mass synchronization for the case N' = 
2 (as well as N' = 3) in Eq. ([ij. In order to explain this 
analytically we rewrite the amplitude equation (jl5p for 
the case TV' = 2 as 

1 - 

h = -lin + i — ^)((Tiiri cos(q;ii) 

+fl(J2ir2 cos(7/'i -i>2+ ai2)), 
1 -\- 

^1 = uji - {— — 3-)(a-iiri sin(aii) 
2ri 

+Aio-2i?'2 sin(V'i - ■02 + "12)), (16) 
1 — 

h = -727-2 + ( ^)(ct22''2 COs(q;22) 

-f (T2l7'i cos(7/;2 - "01 + a2i)), 
1 + 7-"^ 

lp2 = UJ2-{— ^)(cr227'2 Sin(a22) 

2r2 

+o-2i?'i sin(Vj2 - ^1 + "21))- (17) 



77 = 1,2,. ..,A^'. (12) 



Now, substituting the Poisson kernel ([TT|) into the order 
parameter equation © we get 



z„ = / a^{uj'^'^\t)g{uj'-''y)dJ''\ (13) 
' —00 

Assuming a Lorentzian distribution ([5|) for g(w*^''-') and 
evaluating the integral in equation (1131) by contour inte- 
gration we get z* = ari{uj.q — ijri,t)- Thus equation (fT^ 
becomes 



N' 

- E 

=1 



To start with we will consider the case /i = in equa- 
tions (|16l) and (IT7| for analytical simplicity to identify 
pathological synchronization in neuronal populations so 
that the source strongly (completely) drives the target 
populations. However, in general one can observe the oc- 
currence of pathological synchronization for ^ < 1 also 
as in the case of numerical simulations where /x =0.2 and 
0.3. For this purpose, we will also consider below the case 
^ 7^ in equations (ITBl) and (fT7|) so as to provide a theo- 
retical basis for the numerical findings. Consequently, we 
can conclude that when the strength of the synchroniza- 
tion increases in the source it induces synchronization in 
the target. 

From equation (jl6p . for /i = 0, the synchronization of 
the source is characterize d by the stability of the fixed 
point rf = ^1 — 27i/(Tii, where a = cr cos(a). On the 
other hand, the desynchronization state is characterized 
by the stability of the fixed point rf = 0. When cth < 
77 = 1,2,...,A^'. (14) 271, the fixed point rf becomes stable and there is no 
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synchronization in the source. In this state the equation 
for the target population is given as 

1 — 

r2 = -72''2 + ( ^)(ct22?'2 C0s(a22) 

1 + 

V)2 = u)2-{— — ^)cr22r2 sin(Q;22)- (18) 
2r2 

Again one can check that for CT22 < 272 the fixed point 
r2 = is stable and the target population is desynchro- 
nized. Thus when the source is desynchronized, the tar- 
get population is also desynchronized. This result holds 
good for /i < 1 also as is evident from our numerical find- 
ings as depicted in panels (a) and (b) of Figs. [2] and [3] 
for ^ = 0.3. 

On the other hand, when the coupling strength in the 
source an increases so that cth > 271 the fixed point rf 
becomes unstable and rf becomes stable thus establish- 
ing synchronization in the source. The synchronization 
strength of the source increases as yjl — 271/0-11. Af- 
ter synchronization in the source is established, equation 
([T7| reduces to 

1 - 

r2 = -72''2 + ( ^)(0-22?'2 COs(a22) 

+o"2i Vl - 271/CTii cos(Vj + a2i)), 

1 _|_ J.2 

i' = ^ - {—^ — ^)(cr22?'2 sin(Q;22) 
2r2 

+0-21^1 - 27i/o'ii sin(?/' -I- a2i)), (19) 

where ip = ip2 — ipi and ui = UJ2 — UJ1 + (cth —71) tan(aii). 
This equation does not admit the fixed point r2 = 0. 
This means that when synchronization emerges in the 
source, the target population also begins to get synchro- 
nized. The strength of the synchronization in the target 
increases according to 1721 \/l — 271/0-11, eventually lead- 
in g to synchronization with the source (we note here that 
in [3a, |33| an equation similar to ((T^ has been analyzed 
in connection with the stability of the synchronized and 
desynchronized states). This results in the pathological 
strong synchronization which corresponds to our numer- 
ical findings as depicted in panels (c) and (d) of Figs. [2] 
and [3] 

Even for n 0, from equations ((T5)) and (HZ]), we can 
establish that when the source is completely desynchro- 
nized the target is also completely desynchronized. That 
is when the fixed point ri = exists, r2 = also exists. 
When Oil increases, synchronization is induced in the 
source. Hence, the fixed point ri = does not exist as 
a result of which r2 = also becomes nonexistent. This 
means that when synchronization arises in the source, the 
target also spontaneously becomes synchronized. This 
can be described as follows: The eigenvalue of the fixed 
point (ri, r2) = (0, 0) is 

K 1 1 

A± = -7+-±-(^-iylAw-Aa;2)5 -id;, (20) 

where 71,2 = 7^ = 0, i,j = 1,2, k = on + 022, 
A = an- 022, £, = (jA^ + /iofi). Aw = cji - ^2 and w = 



{u!i + uj2)/2. The stability of the fixed point (ri,r2) = 
(0,0) is determined by the real part, i?e(A±) = —7 + 
f ± i(C)^ for Aid = 0, of equation ^ (that is the fixed 
points are stable when Re{X±) < 0). In this connection 
we wish to note that when an = CT22 and fj, — I the 
stabilit y eq uation reduces to the one studied by Montbrio 
et al. j38| for a system of two interacting populations. 
In this case the stability diagram for equation (|20l) will 
be similar to Fig. 1 of 38]. When on increases, the 
condition 7 < ^ ± i(^)2 is satisfied and hence the fixed 
points (ri, r2) = (0, 0) become unstable. This means that 
the increase in an to a sufficiently high value induces 
synchronization in the target also apart from inducing 
synchronization in the source. 

B. Controlling mass synchronization 

In Sec. 3 we have shown from numerical analysis that 
pathological synchronization can be controlled by pro- 
viding a delay feedback. Therefore including the delay 
feedback term in Eq. ([B]) provides the following ampli- 
tude equation. 




where z^ir — Zrji{t — r). In the pathological synchroniza- 
tion state all the populations are in synchronization with 
the source and hence the dynamics of all the populations 
are similar to that of the source and are completely syn- 
chronized with it. Thus all the populations have the same 
synchronization strength and phase, that is ~ i? and 
~ (/), 77 = 1,2,...,A^'. With this reasoning equation 
(PT|) becomes 

1 — i?^ 

R = -7-R+( — ^ — ){iyR~ FRrCos{(l)~(j)r)), 

^ = uj- (i±^)(i>i? - FRr sin((/. - 00), (22) 

for 7 = 7,,, w = w,,, 1/ = S^^i cTw' cos(a,,,,') and 

^ = I],^=i o"!)!)' sin(a,p,')> "H = l,2,...,iV'. Here R 
and 4> are the global order parameters defined by Z = 

Re^"^ = Y^^'=i rjj'e'^'^'^' . Now with a chage of variable 

(j> = ilt + (f) equation (l22|l can be rewritten as 

R = -jR+i — - — )i9R- FRrCOs{(f>-(f>r + nT)), 

4> = {j-{l±^){vR-FR,siu{4)-^r + ^T)), (23) 
Zti 

where w = w — fi. The fixed point i? = of the above 
equation, which characterizes the fully desynchronized 
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state, is not of physical interest here. So we are now 
interested in the pathological synchronization state of the 
system which is represented by the fixed point solution 
i? 7^ and (j) — constant. Consequently from equation 
we have 



y — Fcos(riT) ' 

„ / „ /„ i) — i^sinffir) 

n = oj+{j-9 + Fcos{nT)^— 24 

V — t COS(siT) 

For suitable values of F and r the strength of the patho- 
logical synchronization R can be reduced for a given set 
of system parameters. This is evident from Fig. [5] where 
we have plotted the order parameter i? as a function of 
the strength of the feedback F for two different values 
of the time delay parameter r. From Fig. 8, we observe 
that in the absence of feedback, that is F = 0, there is a 
strong synchronization in the system and hence R ~ 0.9. 
Upon introducing the feedback, the synchronization in 
the system is destroyed for a suitable values of F for a 
given T. The solid and the dotted curves are the analyt- 
ically obtained solutions given by equation (|24l) and the 
curves denoted by * and * are the corresponding numer- 
ical simulations of equation ^ ioi N' = 2. Thus our 
analytical findings are in reasonable agreement with the 
numerical results. 

Thus we find that the above analytical consideration 
correlates our earlier numerical study and provides a 
strong theoretical basis to identify the hypothesis that 
when strong synchronization occurs in the source it in- 
duces synchronization in the target also by increasing the 
coupling strength between the source and the target. We 
further conclude that strong synchronization can be con- 
trolled by a suitable delayed feedback control. 

VI. CONCLUSION 

Occurrence of strong pathological synchronization in 
neuronal populations severely impairs brain function and 
causes pathological tremors. In general, during patho- 
logical synchronization, a particular group of neurons in 
certain region of the brain gets synchronized and acts as 
a pacemaker which drives other cortical neuronal popu- 
lations onto mass synchronization. The synchronization 
frequency in such cases is close to or equal to the tremor 
frequency [l^. 

In this paper we have considered a model of coupled 
populations of phase oscillators to study the occurrence 
of strong synchronization. We consider one of the popu- 
lations to be the source and the other populations that 
are driven to be the targets. We find the occurrence of 
pathological synchronization in the model, and from the 
dynamical properties we are able to hypothesize the pos- 
sible cause for the occurrence of such strong synchroniza- 



tion: When strong synchronization occurs in the source 
this in turn increases the coupling strength between the 
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FIG. 8: The occurrence of desynchronization denoted by a 
decrease in R due to delayed feedback. The solid and the 
dotted curves are the analytically obtained solutions given by 
equation (|24|l . * and ★ curves are the corresponding numerical 
simulations of equation ((!]) for N' = 2, j = 0.02, iu = 1.5 
an — V and Oiij = 0, i, j = 1, 2. 



source and the target. In this state the source strongly 
drives the target and as a result the target gets synchro- 
nized with the source. Thus not only the coupling be- 
tween the source and the target causes the pathological 
synchronization but also the strength of the synchroniza- 
tion in the source plays a vital role. 

We have proposed a possible method to control the 
occurrence of this pathological synchronization. We find 
that a delayed feedback given to the system of neuronal 
populations destroys the strong synchronization and es- 
tablishes desynchronization. The method is demand- 
controlled since the strength and the time delay of the 
feedback is determined from the strength and the fre- 
quency of the synchronization. We have analytically an- 
alyzed the model in the thermodynamic limit and have 
deduced the amplitude equation using the ansatz of Ott 
and Antonsen [361] . From this equation we have shown 
the existence of synchronized solutions in the target when 
strong synchronization occurs in the source. Upon apply- 
ing delayed feedback we have found that the strength of 
the synchronization can be reduced. We believe that our 
hypothesis for the cause of the occurrence of pathological 
synchronization will be useful for a better understand- 
ing and further research in this direction. The proposed 
method will help in the effective stimulation of neuronal 
populations in deep brain stimulation therapy. 

The work is supported by the Department of Science 
and Technology (DST)-Ramanna program (ML), DST- 
IRHPA research project (ML and VKC) and an INSA 
Senior Scientist fellowship (ML). JHS is supported by a 
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